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ABSTRACT 

We present results from multiwavelength observations of the galaxy NGC5005. We 
use new neutral hydrogen (HI) observations from the Very Large Array to examine 
the neutral gas morphology and kinematics. We find an HI disc with a well-behaved 
flat rotation curve in the radial range 20^^140". Ionized gas observations from the 
SparsePak integral held unit on the WIYN 3.5m telescope provide kinematics for the 
central ID”. We use both the SparsePak and HI velocity helds to derive a rotation 
curve for NGC5005. Deep 3.6/rm observations from the Spitzer Space Telescope probe 
the faint extended stellar population of NGG5005. The images reveal a large stellar disk 
with a high surface brightness component that transitions to a low surface brightness 
component at a radius nearly 1.6 times farther than the extent of the gas disk detected 
in HI. The 3.6^m image is also decomposed into bulge and disk components to account 
for the stellar light distribution. Optical broadband B and R and narrowband Ha from 
the WIYN 0.9m telescope complement the 3.6/im data by providing information about 
the dominant stellar population and current star formation activity. The neutral and 
ionized gas rotation curve is used along with the stellar bulge and disk light prohles 
to decompose the mass distributions in NGG5005 and determine a dark matter halo 
model. The maximum stellar disk contribution to the total rotation curve is only about 
70%, suggesting that dark matter makes a signihcant contribution to the dynamics at 
all radii. 

Subject headings: galaxies: kinematics and dynamics; galaxies: structure 
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INTRODUCTION 


Understanding the role and significance of dark matter in the evolution of baryonic components 
(i.e., conversion of the gaseous disk into stars) is a critical aspect for realistic models of galaxy evo¬ 
lution. Observational studies aimed at addressing the influence of dark matter in galaxy formation 
and evolution often use mass decompostion to constrain the distribution of dark matter. Many 
of these projects utilize neutral hydrogen (HI) synthesis observations to derive the kinematics for 
rotation curve decompostion (i.e., The HI Nearby Gal axy Survey (THINGSl: IWalter et al.l 120081 . 


Westerbork Survey of HI in Spiral Galaxies (WHISP); Ivan der Hulst et al.l 120011 ). In this paper, 
we take advantage of the well-behaved outer disk kinematics in the galaxy NGG5005 to perform a 
mass decomposition analysis. We use radio synthesis observations of HI and GO to probe both the 
gas surface density and dynamics as the gas follows the gravitational potential of the galaxy. The 
stellar contribution to the overall mass budget of the galaxy is determined from deep near-infrared 
(NIR) observations. Finally, we use optical broadband B and R and narrowband Ha observations 
to supplement the NIR and radio synthesis data to provide information about the dominant stellar 
populations and star formation activity. 

Large projects like THINGS and WHISP have provided a foundation for numerous other studies 


van Evmeren et al. 

2011. 

Holwerda et al. 

2011, 


Vollmer &: Lerovll201lh and dark matter haloes fe.g.. lOh et al 


2011 


de Blok et al.ll2008l L However, 


an examination of the statistical significance of connections between the dark matter and baryonic 
properties of these galaxies may be biased due to the representative nature of the samples. In an 
effort to address fundamental questions regarding the growth and distribution of stellar disks in 
dark matter haloes in a statistical manner, we have undertaken a project correlating structural 
properties and star formation activity with the dark matter properties of the host galaxy. The 
project uses a statistical sample of 45 nearby galaxies and builds on ex isting data obtained as part 
of the Extended Disk Galaxy Exploration Science (EDGES) Survey (jvan Zee et al.l 120121 ). which 
includes deep Spitzer Space Telescope imaging observations. The statistical sub-sample consists 
of all galaxies in the EDGES Survey that are optimally suited for rotation curve decomposition, 
based primarily on the limited range of inclination angles that yield both accurate rotation curves 
and surface density profiles. Our combination of optical, NIR, and radio synthesis observations 
will enable a comprehensive statistical analysis of the evolution of galactic disks as a function of 
baryonic mass, environment, and dark matter mass. 

In this paper, we focus our efforts on a single galaxy from the sam ple, NGC5005. This 


SAB(rs)bc type galaxy has been studied extensively in th e literature (e.g., iBurbidge et al 


1961 


Blackman 

1979. 

Duric et al. 

1996. 


Batcheldor et al 


nuclear emission-line region (LINER) in its center (jHo et al.l Il997u and has been classified as a 


200511. It is k nown to host a low-ionization 


low-luminosity active galactic nucleu s (LLAGN) with a polycyclic aromatic hydrocarbon (PAH)- 
dominated mid-infrared spectrum bv iMason et al.l (j2012l L The morpholo gy of NGC5005 as vi ewed 
in the infrared and optical features a strong bulge with boxy isophotes (lEskridge et al.ll2002l ) and 


a complicated dust morphology with no coherent structure (iHughes et al 


2003). There is also 











































































- 3 ~ 


evidence of a fairly weak, rou nd bar whi c h is n early aligned with the major axis of the galaxy. De¬ 
spite the weakness of the bar, iDas et ahl (|2003l i found that it could still be responsible for creating 
shocks visible in the central m olecular gas velocity fie ld. The molecular gas dynamics in NGC5005 
were studied in more detail bv I Sakamoto et ahl (|2000l l. Their observations reveal CO emission from 
the nucleus, a ring at about 3 kpc, and a stream of gas northwest of the nucleus. The S-shape 
in the nuclear velocity field in 


Sakamoto et al 


([20001) is a signature of non-circular gas motions. 
Sakamoto et al.l (120001 ) cite a stellar bar of length ~5 kpc as being accountable for infalling gas 


seen as the northwest stream and for creating the 3 kpc ring which may form at the ultraharmonic 
resonance of the bar or alternatively constitute a pair of spiral arms that originate at the bar’s 
ends. A position-velocity cut along the major axis of NGC5005 reveals a steep rise in the rotational 
velocity of the CO emission (their Figure 5). The rotational velocity rises to approximately 300 
km s“^ within the central lO'h 


NGC5005 lies in a relatively populated environment that is near the Ursa Major Cluster. The 
distance to NGC5005 is uncertain; distance measurements rang e from 13.7 Mpc to 34 .6 Mpc. In 
this paper, we adopt the Type la supernova distance of 16.5 Mpc ([Takanashi et al.ll2008l). NGC5005 


is tho ught to be part of a physical galaxy pair with the nearby spiral galaxy NGC5033 (jHelou et al 


1982l l. The two galaxies are separated by a projected distance of approximately 41' (~200 kpc at 


16.5 Mpc) and are at similar redshifts. There are no clear signs of any interaction between the galaxy 
pair. There are two dwarf galaxies, SDSS J131115.77-1-365911.4 and SDSS J131051.05-1-365623.4, 
that are within 10' of NGC5005 and also lie at similar redshifts. Farther out, there are many 
additional dwarf galaxies within one degree of NGC5005 which are part of the NGC5033 group. 


The goal of the present work is to investigate the dark matter and baryons in NGC5005. 
NGC5005 is also a proof of concept for the larger study. The observational data are presented 
in Section [2] and we focus on the kinematics in Section [3j Section U] describes the rotation curve 
decomposition analysis and the determination of the dark matter distribution in NGC5005. Section 
[5] provides a brief summary of the results. 


2. OBSERVATIONAL DATA 

Multiwavelength observations have been acquired to probe the stellar and gas content in 
NGC5005. Deep NIR images taken at 3.6//m from the Spitzer Space Telescope trace the extended 
stellar population, while moderate depth optical broadband B and R and narrowband Ha provide 
information about the dominant stellar population and star formation activity. Radio synthesis 
observations of neutral hydrogen from the Very Large Array (VLA) 0 are used to probe both 
the gas surface density and kinematics. Spectroscopic integral field unit (IFU) observations from 


^The Very Large Array is operated by the National Radio Astronomy Observatory, which is a facility of the 
National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
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SparsePak on the WIYN 3.5m telescope provide ionized gas kinematics in the central region of 
NGC5005. Figured] shows the new VLA HI observations along with optical and NIR imaging data. 
Table dl provides a summary of integrated properties derived from this multifrequency dataset. We 
describe the data acquisition and processing below. 


2.1. VLA Observations 

Radio synthesis observations of NGC5005 using the 21 cm line of HI were obtained with the 
VLA in C configuration on 2013 August 4 and 5 for a total on-source time of 407 minutes. The 
data had an initial frequency resolution of 7.8125 kHz channel”^ (~1.65 km s“^ channeH^). The 
standard flux calibrator 3C286 was observed at the beginning and end of each observing block, and 
the phase calibrator J1227-I-3635 was observed approximately every 40 minutes, so that the data 
may be flux and phase calibrated. 

In order to flag uniformly, the data were loaded into AIPs[^ and the inner 75% of the bandpass 
was combined to generate a “channel zero” dataset. The channel zero data was flagged and phase 
calibrated. The calibration solutions were applied to the line data, which was then bandpass 
calibrated with 3C286. After calibrations, the data were continuum subtracted in the uv plane and 
then combined. The combined, continuum subtracted datasets were imaged using the AIPS task 
IMAGR. Three cubes were created with different Robust weighting parameters for varying spatial 
resolution and a channel averaging of 3 for a resulting velocity resolution of ~5 km s“^ channel”^ 
(for simplicity, these cubes are designated as low, medium, and high; Table [2|). Only the cubes 
with the highest sensitivity (low) and highest resolution without severe loss of sensitivity (high) are 
shown in Figured] A lower velocity resolution cube (15 km s“^ channel”^) was additionally created 
to probe the lowest column density of HI gas. This velocity resolution was chosen to match the 
channel width to the expected line width from thermal gas motions. This cube was also smoothed 
spatially post-imaging to increase the sensitivity to low column density gas. 

In addition to the new VLA C configuration data, there are existing archival VLA D config¬ 
uration data taken on 2000 August 11 for a total on-source time of 80.3 minutes (project code 
AW536). These archival data were retrieved from the archive and processed in AIPS. Properties 
of the imaging of this data set are presented in Table [2] The archival D configuration data was not 
combined with the new G configuration data due to the lower velocity resolution and the smaller 
frequency coverage of the archival data. The galaxy signal nearly fills the bandpass in the D config¬ 
uration observations. Despite this, the integrated HI flux and HI mass measured from the archival 
data is consistent with the new C configuration values within the measurement errors. Both this D 
configuration data and the spatially smoothed, 15 km s“^ channel”^ velocity resolution cube were 
used to confirm the lack of detection of low column density gas in the higher resolution cubes. 


^The Astronomical Image Processing System (AIPS) has been developed by the NRAO. 
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Fig. 1.— (a) Low spatial resolu tion HI integrated intensity contours (white) and CO intensity 
contours from lHelfer et al.l (j2003l l (red) overlaid on the Spitzer 3.6p,m image. Note that this field of 
view focused in on NGC5005 does not include the entire low surface brightness stellar disk, (b) Low 
spatial resolution HI intensity contours overlaid on the low spatial resolution HI intensity image. 
Note the hole in HI emission in the center, (c) High spatial resolution HI intensity contours overlaid 
on the Ha image, (d) HI velocity field, or first moment map with isovelocity contours derived from 
the low spatial resolution HI data. The velocity contour values are indicated in the plot in units 
of km s“^. The HI intensity contours represent column densities of 10^® atoms cm“^ x (1, 2, 4, 8, 
12, 16) in images (a)-(c). 
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Table 1. Integrated Properties of NGC5005 


Property 

NGC5005 

Distance®’ 

16.5 Mpc 

Morphological Type'^ 

SAB(rs)bc 

UlB 

10.74±0.02 

mR 

9.29±0.02 

ms.e 

6 . 21 ± 0.01 

Mb 

-20.50±0.02 

(B-R)o 

1.43±0.02 

(B-3.6)o 

4.48±0.02 

(R-3.6)o 

3.05±0.02 

log(Ha fiux) 

-11.42±0.09 erg cm“^ 

Equivalent width 

7.93±0.32 A 

Total SFR 

0.67±0.14 M© year-1 

Total stellar mass*^ 

(9.16±1.83)xl0i“ Mq 

HI flux 

19.0±3.8 Jy km s"! 

Total HI mass 

(1.22±0.24)xl0® Mq 

CO flux^' 

1278±484 Jy km s"! 

Total H 2 mass‘d 

(2.73±1.03)xl0® M 0 

Total dynamical mass® 

(2.82±0.47)xl0ii M 0 


Note. — The apparent magnitudes are measured values and are not cor¬ 
rected for extinction. The reported apparent magnitudes and colors are mea¬ 
sured at R 25 . Mr and colors are extinction corrected assuming Ab=0.051 
and Ab= 0.031 llSchlaflv & Finkbeiner|[201lh . The extinction correction for 
the NIR is assumed to be negligible. The Mr value also includes a small 
correction ( 0.10 mag) based on an extrapolation of the observed stellar disk. 

^Takanashi et all ll 200 g|) 

^RC3; Ide Vaucouleurs et al. 1 119911) 

^Assumed best-fitting M/L 3.6 = O.SdiO.l 

* 1 Heifer et al.l 1 I 2 OO 3 II : mass corrected for distance=16.5 Mpc 

^Measured at Rbreak? assuming the rotational velocity stays constant. 


Table 2. HI Synthesis Image Parameters 


Image Name 

Velocity Resolution 
(km s“^) 

Robust 

Beam Size 
(arcsec) 

Beam Position Angle 
(degrees) 

Noise 

(mjy beam“^) 

low 

5.0 

5 

26.4x17.2 

64.2 

0.57 

medium 

5.0 

0.5 

20.5x14.8 

67.0 

0.61 

high 

5.0 

- 0.5 

16.9x13.2 

72.6 

0.72 

smooth 

14.9 

5 

52.9x34.5 

64.2 

0.53 

D config^ 

10.4 

5 

59.6x51.8 

-45.2 

0.96 


Archival VLA D configuration (project code AW536) 
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The total intensity map and velocity field of the final combined C configuration data imaged 
at the low spatial and 5 km s“^ channel”^ velocity resolution are presented in Figures [T] (a), (b), 
and (d). The velocity field will be discussed in greater detail in Section [3Tl Figure [I](c) shows the 
same data in contours imaged at high spatial resolution. The first contour in Figures [T] (a), (b), 
and (c) corresponds to a column density of 10^^ atoms cm~^. There are two notable features of the 
HI emission in NGC5005 visible in Figures [T] (a) and (b). First, there is a slight hole or depression 
in the HI column density at the center. This is expected to occur in galaxies where there is a phase 
change from atomic to molecular gas in the central kilo parsec. Bar-driven transp ort also tends to 


conce ntrate the molecular gas in the centers of galaxies. (|Sheth et al.ll2005l : see also lSakamoto et al 


2000 1. Archival CO data allows us to trace the molecular gas component in the nucleus to fill in 


the missing gas surface density information (Section l2.2p . Second, the gaseous disk does not extend 
nearly as far as the stellar disk. Indeed, NGC5005 is kno wn as an HI-deficient galaxy with a 
measured mass in HI that is almost half of its expected value (jFumagalli et al.ll2009l l. Furthermore, 
we did not find low column density gas in either the lower velocity resolution, spatially smoothed 
map or in the archival D configuration data beyond what is shown in Figured] (a). The higher 
resolution HI contours displayed in Figured) (c) show that the neutral gas appears to be spatially 
coincident with the ionized gas seen in Ha, as is expected for star forming regions in galaxies. 
Global properties of NGC5005 measured from this HI dataset are presented in Tabled! 


The field of view of the HI observations includes not only NGC5005, but also nearby dwarf 
galaxies SDSS J131051.05+365623.4, SDSS J131115.77+365911.4, SDSS J131105.57+371036.1, and 
SDSS J131126.81+371842.3. Only SDSS J131051.05+365623.4 is detected in the HI images (Figure 
[2]) with a measured flux of 0.21 Jy km s“^ and line width at 50% of the peaks of 33 km s“^; this 
yields a corresponding HI mass of 1.33x 10^ Mq at the adopted distance of 16.5 Mpc. Outside of the 
primary beam, we also detect NGC5002 and NGC5033. No other known or suspected companions 
were detected in the HI observations. 


2.2. Archival CO 


The central hole in HI column density makes it difficult to extract trustworthy total gas 
surface density measurements from the VLA data alone. Therefore , archival CO data from the 
BIMA Survey of Nearby Galaxies (BIMA SONG) ( Heifer et al. 2003 ) was acquired to supplement 
the HI observations of NGC5005. Observations for the BIMA SONG survey were carried out from 
1997 November through 1999 December using the Berkeley-Illinois-Maryland Association (BIMA) 
millimeter interferometer in Hat Creek, CA. The archival CO data cube for NGC5005 has a velocity 
resolution of about 10 km s“^ and a spatial resolution of 6.2'' x 6.0". The integrated CO intensity 


i s ma pped in red contours in Figure [D (a). A more detailed map can be found in iHelfer et al 


(120031 1. This map allows the central molecular gas surface density to be determined. Note that 


the spatial extent of the molecular gas is confined within the central ~40" or ~3.2 kpc at 16.5 
Mpc. As is typical for barred galaxies with CO emission, the molecular gas is concentrated in the 
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Fig. 2.— Low spatial resolution HI integrated intensity contours (white) overlaid on the Spitzer 
3.6/im image to show the full extent of the stellar disk. The 3.6/im image on the left is smoothed 
and edited to highlight the low surface brightness features, including a possible dwarf companion 
to the north and a hint of a potential stellar stream connecting it to NGC5005. The outlined 
region demarcates where flux assumed to be part of the potential stream was measured. The image 
on the right shows the unmasked, unsmoothed large field of view around NGC5005. The dwarf 
companion SDSS J131051.05+365623.4 (labeled as SDSS 1310+3656 in the left panel) to the south 
of NGC5005 is detected in HI. The dwarf companion SDSS J131115.77+365911.4 (abbreviated as 
SDSS 1311+3659) is also visible to the southeast of NGG5005 in the left panel. The white dashed 
ellipse indicates Rbreak (Section [2+]) in both images. 
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center (jSheth et alJl2005l ). In the present study, we focus on using the archival CO observations 
to improve the accuracy of the gas surface density measureme nts only. A complete a nalysis of the 
complex molecular gas dynamics in NGC5005 can be found in ISakamoto et alJ (2000). 


2.3. WIYN SparsePak 


We obtained an H a velocity field for the central region of NGC5005 using the SparsePak IFU 
(jBershadv et al.ll2004l ') on the WIYN 3.5m telescope. The SparsePak IFU is composed of eighty-two 
f)'' diameter fibers arranged in a fixed 70" x 70" rectangle. Observations were acquired on the night 
of 26 April 2014 during reasonably clear, but non-photometric weather with good seeing. We used 
the STAl CCD binned 4x3 with the 316@63.4 grating in eighth order, centered at 6686 A to 
cover a wavelength range from 6480 A to 6890 A. This resulted in a velocity resolution of 13.9 km 
s“^ pixel”^, which is sufficient to derive an accurate velocity field for NGC5005. 

For simplicity, the SparsePak array was aligned on the sky with a position angle of 0°. We used 
a three pointing dither pattern to spatially fill in gaps between fibers. We took three exposures at 
900 s for each dither pointing in order to be able to detect diffuse ionized gas, not just star forming 
regions. Additional observations of blank sky were also taken to remove sky line contamination 
more accurately, as NGC5005 is much more extended than the SparsePak held-of-view (see Figure 

ED. 

The SparsePak data was processed using standard tasks in the HYDRA package within irafI 
The data were bias-subtracted and flattened, and the IRAF task DOHYDRA was used to fit and 
extract apertures from the IFU data. The spectra were wavelength calibrated using a wavelength 
solution created from ThAr lamp observations. The three exposures for each dither pointing were 
cleaned of cosmic rays. The individual images were sky subtracted using the separate sky pointing 
scaled to the strength of the 6577 A sky line, which is close to the redshifted Ha and [Nil] A6584 
lines in some locations. The cleaned, sky subtracted images were averaged together to increase the 
signa l-to-noise r atio and then flux calibrated using observations of spectrophotometric standards 
from Oke (1990). 


We illustrate representative spectra from four SparsePak hbers in the central region of NGC5005 
in FigureElto demonstrate the types of spectral emission found in NGC5005. The locations of these 
fibers start in the central nucleus of NGG5005 and move out towards the northwest just off the 
minor axis (see the rotated Ha image on the right). While our primary analysis relies on cross¬ 
correlation of the entire spectrum, we also fit Gaussian profiles to emission lines of Ha, [Nil], and 
[SH] in all the spectra to estima te the full width half maxim um (FWHM) values and the emission 
line ratios. As seen previously (jBurbidge &: Burbidgd 119651 1 the strength of the [Nil] A6584 line 


®IRAF is distributed by NOAO, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 
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emission relative to Hq in the central fiber (+000+000) implies non-thermal emission processes 
([Baldwin et alJll98lh . The broadness of the nuclear e mission lines fr om the center of NGC5005 
are suggestive that the nucleus is accretion powered (|Ho et alJ 1199711 . We find an approximate 
FWHM of 900+200 km s“^ based on single Gaussian fits to [Nil] and [SII] lines within the cen¬ 
tral 5" fiber. The second spectrum in Figure [3] (-005+003) was extracted from the adjacent fiber 
to the northwest of the central fiber. This spectrum displays double-peaked em ission lines that 


sugge sts either infalling or outflowing gas. Indeed, this fiber is placed near where I Sakamoto et al 


(|2000l l detected infalling molecular gas in their CO observations of NGC5005. The third fiber’s 
spectrum (-010+006) still shows broad emission lines with non-thermal line strengths, but less so in 
both aspects than the central fiber. This spectra is representative of diffuse ionized gas emission in 
NGC5005, which we measure to have typical FWHMs of about 140+10 km s“^ and [Nil] A6584 /Hq! 
line ratios of 3.0+1.2, which is consistent with non-thermal or shock excited gas. The fourth fiber’s 
spectrum (-015+008) shown in the bottom panel of Figure [3] is placed on the ring of Ha emission 
surrounding the nucleus. This spectrum displays narrower line widths and stronger Ha to [Nil] 
A6584 line emission strengths. We find typical FWHMs of approximately 80+5 km s“^ and [Nil] 
A6584/Ha ratios of 0.4+0.1 in similar spectra, which are ii idicative of primarily thermal emission 
processes most likely due to recent star formation activity ([Baldwin et al.lll98ll ). 


Inspection of spectra extracted from a few other fibers revealed additional locations with 
double-peaked emission lines, similar to the -005+003 fiber spe ctrum. These sp e ctra w ere located 
approximately east-west on the sky along where evidence from I Sakamoto et al.l ([2000l ) suggests a 
stellar bar exists. We note, however, that double-peaked emission lines may also occur due to 
smearing of the rotation curve within the f)" SparsePak fibers. The spectra along the major axis 
would be most susceptible to this beam smearing as the velocity gradient is highest here. Since 
the bar is oriented close to the major axis of NGG5005, it would be difficult to verify that the 
double-peaked emission lines are signatures of infalling gas along the bar. Alternatively, the beam 
smearing effect would be minimized for velocities along the minor axis, so it is reasonable to assume 
that the double-peaked emission lines seen in the -005+003 spectrum in Figure [3| are real and are 
not merely due to observational effects. 

The flux calibrated, sky subtracted spectra were cross-correlated with a template emission line 
spectrum to extract the luminosity weighted mean recessional velocity at every position using the 
IRAF task FXGOR. We used the spectrum of an HII region from within NGG5005 as our reference 
template for the cross-correlatio n analysis due to its strong narrow emission lines. Spectra with 
a Tonry and Davis ratio (TDR; iTonrv &: DavisI 119791 ) less than 5 were not included in deriving 
the velocity field. Spectra with broad or double-peaked emission line profiles resulted in broad 
cross-correlation peaks that were nonetheless of high significance and had a well-defined centroid. 
The velocity values from FXGOR were then placed into a grid that mapped the SparsePak fiber 
locations. To fill in the missing spacings, the velocity field was then interpolated using the average 
value from the nearest eight pixels. Analysis of the velocity field will be discussed further in Section 
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Fig. 3.— Example spectra extracted from four different fibers from the SparsePak IFU. The 
positions of the SparsePak fiber apertures on NGC5005 are displayed on the narrowband Ha image 
from HDI on the right. The fibers from which the four spectra were extracted are highlighted in 
white. The image has been rotated such that north is to the right and east is up. The numbers in 
the top left corner of each spectrum panel indicate the fiber’s position offset from the central fiber 
in arcseconds. 
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2.4. Spitzer 3.6/im Data 


Imaging observations at 3.6//m were obtained using the Infrared Array Camera (IRAC) on 
the Spitzer Space Telescope on 2012 July 12 and 19. Observations were conducted in a grid-like 
mapping pattern of 10 x 9 dithered pointings, where each pointing was offset by ~100'' to cover 
a total area of 18' x 33'. Observations were taken in two sets separated by several days to enable 
asteroid removal from the mosaicked images. The resulting map has a total exposure time of 1800 
s pixel“^. This exposure time provides a sensitivity of a fewx 0.01 Mq pc“^ over a held of view 5 
times the size of the optical radius R 25 . After standard processing, the individual pointings were 
combined into maps using the MOsaicker and Point source Extractor (MO PEX) software. Images 
were drizzled to resample the data to 0.75" pixel”^ ( Fruchter &: Hookll2002l l. The maps were ht for 
a hrst order sky subtraction to remove any gradient in the sky level. 

Figure [2] shows the full extent of the faint stellar disk component of NGC5005 detected in the 
hnal mosaicked and sky subtracted 3.6/im image. The image on the left has been smoothed to 
bring out the lowest surface brightness features. The immense sensitivity over such a large held 
of view enables a complete study of the faint extended stellar component of NGC5005. NGC5005 
has a large high surface brightness disk, which is surrounded by an even more extended low surface 
brightness component. This faint stellar population extends about twice as far as the break or 
transition radius between the high and low surface brightness components and nearly three times 
farther than the gas disk detected in HI. This is an unusual result, as most spir al and irregular 
galax ies with HI have a diameter in HI that is larger than the optical diameter (|Broeils &: Rhee 

I997I 1. 


There is also a feature that appears to be a faint stellar stream extending to the northeast 
towards an apparent companion galaxy SDSS J131105.57-1-371036.1. To measure the stellar mass 
associated with the apparent stellar stream, we used the smoothed and masked image to create a 
hand-drawn polygonal aperture that traces the area associated with the diffuse light (Figure [2]). 
We measured a lower limit for the flux in the stream of 0.15 ± 0.13 mJy based on a conservatively 
masked image. However, this value significantly underestimates the total flux associated with the 
stream due to the significant masking required to remove the large area associated with the bright 
foreground star located in the middle of the stream. Thus, to obtain a more accurate estimate of the 
total mass associated with the stream, we replaced the masked pixels with a value interpolated from 
neighboring pixels and remeasured the total flux to be 0.35 mJy. Assuming the best-fitting disk 
mass-to-light ratio of 0.5 (Section H]) and a distance to NGC5005 of 16.5 Mpc, this corresponds to a 
total stellar mass for the stream of 1.8 x 10 ^ Mq, which is an order of magnitude l ess massive than 


notable streams iden tified around other galaxies, such as M63 (iGhonis et al. 


2011 


Staudaher et al 


in prep.h and M83 (jBarnes et al.ll2014l l. The apparent companion galaxy is not detected in the 
HI observations and does not have a redshift measurement available in the literature. However, if 
SDSS J131105.57-1-371036.1 is at the same distance as NGG5005 and we again assume a mass-to- 
light ratio of 0.5, it would have a total stellar mass of 2.3 x 10^ Mq (Table [3]). Given the slightly 
larger mass of the apparent companion galaxy, interactions between the apparent dwarf companion 
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and NGC5005 cannot be immediately ruled out as an origin for the potential stellar stream. 


The 3.6//m image of NGC5005 was decomposed into bulge and disk components to represent 
the stellar mass distributions in the rotation curve decomposition analysis. Note that the emission 


at 3.6 /im is predominantl y stellar, and not greatly affe cted by PAH emission at 3.3^m (iMeidt et al 


20121 ). We used DiskFit ( Spekkens fc Sellwoodll2nn7l ) to derive a photometric model that fit bulge 
and disk components to the 3.6//m image itself. We did not include a bar component in the 
photometric model, despite evidence for its existence (jPas et al.ll2003l . ISakamoto et al.ll 2000 l b since 
we only fit stellar bulge and disk components in the mass decomposition (Section 0]). Our models 
are not sensitive to the exact inner structure of the galaxy, so the lack of a bar component will have a 
negligible effect on the derived dark matter profile. Radial surface brightness profiles of the resulting 
bulge and disk models are shown in Figure |4] (a). We additionally extracted a surface brightness 
profile from the 3.6/rm image using concentric, con-eccentric ellipses to provide a consistency check. 
The DiskFit model and the ellipse photometry yielded similar total surface brightness profiles. 
NGG5005’s surface brightness profile displays a relatively small bulge component with a complex 
disk component featuring a break, or change in slope of the exponential disk where the transition 
from high to low surface brightness stellar disk occurs. The break radius was estimated by fitting 
two exponentials to the radial ranges of 80''-180" (inner) and 280"-443" (outer). The radial ranges 
were chosen so as to avoid locations where the surface brightness profile displays a change in slope. 
Specifically, the inner limit of 80'' avoids the bulge component and the dip in surface brightness 
around 40". Scale lengths of these inner and outer exponential fits are given in Table 01 The 
break radius is defined as the radius at which the two different exponential fits intersect, which is 
216.3" (or 17.3 kpc at 16.5 Mpc). This break radius is indicated on the bottom axis of the panels 
in Figure m along with the optically defined R 25 (see Section 12.5h and radial extent of the HI disk. 
The size of the neutral gas disk is only about 60% of the size of the high surface brightness stellar 
disk, as measured by Rbreak- Figure 0] (b) shows the radial trends in B-3.6 and R-3.6 colors. The 
colors display nearly identical trends with radius featuring a relatively complex distribution in the 
inner ~80" most likely due to star formation in spiral arms. 

Aperture photometry on the 3.6/rm image was carried out using the GALPHOT0 package for 
IRAF. Concentric, con-eccentric ellipses were used to measure the integrated magnitude at R 25 
(Section 12.5p given in Table [H No foreground extinction correction was used for the NIR since the 
expected value is negligible. The total stellar mass presented in Table [U was determined from the 
total stellar luminosity at 3.6/im assuming a mass-to-light ratio (M/L) of 0.5 (see Section 0|). The 
luminosity was calculated from the absolute magnitude at 3.6/im, which includes a slight correction 
(0.04 mag) based on an extrapolation of the 3.6^m surface brightness. The ellipse photo metry was 
also used to estimate the concentration of light in NGC5005 using the definition from IConselice 


^GALPHOT is a collection of scripts in the IRAF STSDAS environment first developed by W. Freudling and J. 
J. Salzer. The current version has been further enhanced by members of the Cornell Extragalactic Group and is 
maintained by M. P. Haynes. 
















Table 3. Properties of Known and Suspected Companions to NGC5005 


Name 

RA, Dec 
(J2000) 

Systemic^ 
Velocity 
(km s“^) 

Angular 

Separation 

(arcmin) 

Physical^ 

Separation 

(kpc) 

Size'^ 

(axb) 

("x") 

Position 

Angle 

(°) 

mB 

mR 

m3.6 

(B-R)o 


NGC 5005 

13:10:56.2, 

37:03:33 

948 

_ 

_ 

336.4x162.4 

65.55 

10.74+0.02 

9.29+0.02 

6 .21+0.01 

1.43+0.02 

-20.50+0.02 

SDSS J131115.77+365911.4 

13:11:15.7, 

36:59:14 

956 

5.85 

28.1 

42.4x21.6 

-22.24 

18.56+0.07 

17.53+0.07 

15.75+0.01 

1.01+0.09 

-12.58+0.07 

SDSS J131051.05+365623.4 

13:10:51.0, 

36:56:23 

1026 

7.24 

34.7 

43.8x28.2 

32.21 

17.76+0.05 

16.80+0.05 

14.81+0.01 

0.94+0.07 

-13.38+0.07 

SDSS J131105.57+371036.1 

13:11:06.2, 

37:10:40 

- 

7.29 

35.0 

62.8x45.4 

-62.08 

17.71+0.10 

16.75+0.10 

15.15+0.01 

0.94+0.14 

-13.43+0.10 

SDSS J131058.75+364943.8 

13:10:58.6, 

36:49:42 

- 

13.72 

65.9 

59.4x36.2 

-49.14 

18.28+0.12 

17.29+0.12 

- 

0.97+0.17 

-12.86+0.12 

SDSS J131058.04+364812.4 

13:10:58.1, 

36:48:12 

- 

15.35 

73.7 

39.6x25.0 

-8.15 

19.46+0.13 

18.25+0.13 

- 

1.19+0.19 

-11.67+0.13 

SDSS J131126.81+371842.3 

13:11:26.9, 

37:18:42 

960 

16.33 

78.4 

63.0x50.8 

-7.90 

17.21+0.06 

16.11+0.06 

- 

1.08+0.09 

-13.93+0.06 

kkh 081 

13:11:11.6, 

36:40:46 

1032 

22.94 

110.1 

58.6x35.8 

17.85 

18.70+0.16 

17.68+0.16 

- 

1.00+0.24 

-12.44+0.16 

NGC 5002 

13:10:37.9, 

36:38:03 

1091 

25.78 

123.7 

212.8x114.4 

-6.64 

14.16+0.04 

13.33+0.04 

- 

0.81+0.05 

-16.98+0.04 


Note. — The apparent m agnitudes are measured values and are not corrected for extinction. B-R and Mb are extinction corrected assuming Ab= 0.051 and Ar= 0.031 
llSchlaflv Sz Finkbeiner|[201lh . The extinction correction for the NIR is assumed to be negligible. I 

h-i 

^From NED, if known (http;//ned.ipac.caltech.edu/ ) ^ 

^Assuming a distance of 16.5 Mpc 

^ Sizes indicate the diameters in arcseconds of the apertures used to measure the reported magnitudes, except for NGC5005 which is the size measured at the optical diameter 
D25- 
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Fig. 4.— (a) Surface brightness profile at 3.6/rm decomposed into bulge (squares) and disk (stars) 
components. The dashed-dotted and dashed lines are exponential fits to the radial ranges 80"- 
180" and 280"-443", respectively. The arrows represent the radial extent of the HI disk, the optically 
defined R 25 , and the approximate radius at which there is an observed change in slope of the outer 
exponential disk, (b) Radial trend in B-3.6 (squares) and R-3.6 (triangles) color, (c) Radial trend 
in B-R color, (d) Radial trend in the log Ha equivalent width. 
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(|2003l l (see Tabled]). Our calculated value of 4.3 is towards the higher end of c oncentrations an d 
is in agreement with th e value of 4.16 measured in the R-band for NGC5005 bv IConselicel (|2003l l. 
McDonald et al. ( 2009l i performed a study of NIR surface brightness distributions of Virgo cluster 


galaxies and discovered a bimodality of concentration with peaks at C = 3 and 4.5 in high surface 
brightness disk galaxies. The concentration value determined here for NGC5005 is consistent with 
the second peak in the bimodal distribution. 


2.5. Optical Data 

Optical broadband B and R and narrowband Ha images provide insight into changes in the 
age or metallicity of the dominant stellar population in NGG5005. Observations were taken with 
the WIYN 0.9m telescope at Kitt Peak National Observatory in February 2013 with the S2KB 
imager and April 2014 with the Half Degree Imager (HDI). NGG5002, which was detected outside 
the primary beam in the VLA data, was also observed with HDI. During both observing runs, 
narrowband imaging was done using a filter with FWHM of about 60 A centered at a wavelength of 
6580 A. A filter with similar FWHM and slightly offset central wavelength was used for continuum 
subtraction. The Ha images were taken with total exposure times of 2 x 20 minutes in the 6580 
A filter, and 20 minutes in the narrowband continuum filter. Broadband imaging was done with 
three exposures of 900 and 300 s with the B and R filters, respectively. An additional short R-band 
exposure was acquired to replace saturated pixels in the core of NGC5005. 

The optical images were reduced and analyzed with the IRAF package. The image reduction 
included bias subtraction and flat fielding. Both the S2KB and HDI images were flat-fielded using 
twilight flats for the broadband filters and dome flats for the narrowband. Prior to combination, 
sky values were measured as the mode of galaxy-free regions of the images. The images were then 
scaled, aligned, and averaged together to make a final combined image. The WCS was updated 
using star lists from the USNO catalog. The S2KB and HDI images were processed independently 
in order to provide internal consistency checks for the derived structural parameters. Since the HDI 

Table 4. Radial Properties of NGC5005 


Parameter 


Value 


Inner scale length (3.6/rm) 
Outer scale length (3.6/im) 
Outer break radius (3.6/im) 
D25 

HI size at 10^° atoms cm“^ 
Concentration (C=51og 
log (Equivalent width) gradient 
Color gradient (B-R) 


31.6” (2.5 kpc at 16.5 Mpc) 
124.6” (10.0 kpc at 16.5 Mpc) 
216.3” (17.3 kpc at 16.5 Mpc) 
336.4” (26.9 kpc at 16.5 Mpc) 
276” (22.1 kpc at 16.5 Mpc) 
4.3 

-0.11 log (EW) arcmin”^ 
-0.13 color arcmin”^ 
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images have slightly higher quality and a larger field of view than the S2KB images, the subsequent 
analysis will focus on results from the HDI observations. 


In addition to NGC5005, our HDI observations include eight known or suspected companion 
galaxies, including NGC5002 (Table [3]). Velocities presented for the known companions in Table 
[3] are from the NASA/IPAC Extragalactic Database (NED)I§ The three unconfirmed companions 
are diffuse, irregular galaxies that do not appear to be larger background galaxies based on visual 
inspection of the broadband images. These galaxies are identified using the Sloan Digital Sky 
Survey (SDSS) designation that most closely matched the photometric center of the galaxies. No 
corresponding HI or Ha emission was detected for the unconfirmed companion galaxies that look 
like they might be associated with NGG5005. The absolute B-band magnitudes presented in Table 
[3] are calculated assuming a distance of 16.5 Mpc. 


Aperture photometry of NGG5005 using concentric, con-eccentric ellipses was carried out to 
derive surface brightness magnitudes for both the B- and R-band HDI images. Table [T] gives 
measured integrated magnitudes derived from the aperture photometry. The integrated apparent 
magnitudes tabulated here are measured at the radius at which the B-band surface brightness 
equals 25 mag arcsec“^ (R 25 ) and are not cor rected for foreground extinct ion. The absolute B-band 
magnitude has been corrected for extinction (jSchlaflv fc Finkbeineiil2nill ) and includes a correction 
(0.10 mag) based on an extrapolation of the B-band surface brightness profile. Integrated colors 
are also measured at R 25 and have been extinction corrected. 


Inspection of the Ha image indicates that there are no visible star formation regions beyond 
R 25 . Aperture photometry was performed on the Ha image to measure an Ha flux and derive 
an equivalent width. The current total star for mation rate (SFR) i s deri ved from the Ha flux 
(measured at R 25 ) using the calibration given in iKennicutt &: EvansI (j2012l l. The calculated SFR 
of 0.67=fc0.14 Mm vr~^ is much lower than the median value of about 1.5 Mq yr“^ found by 


■Tames et al.l (120041 ) for galaxies of the same Hubble T-type as NGG5005. However, this SFR is 


most likely an underestimate due to internal extinction, which we cannot accurately account for 
without a spatially resolved total infrared map due to the strong active galactic nucleus (AGN) 
contamination. The equivalent width (EW) of the Ha emission line is used as a tracer of the 
specific star formation rate. It is calculated by dividing the Ha flux by the continuum flux density 
measured from the Ha narrowband continuum filter. It, therefore, serves as an indicator of the 
strength of the current SFR relative to the past average SFR. A larger EW value would indicate 
a larger current SFR relative to the continuum, or past average star formation. NGG5005 has a 
relatively small EW value of 7.93±0.32 A, which is consistent with the little Ha emission seen in 
the narrowband image. This EW is also more of an upper limit as it is contaminated by strong 
[Nil] emission from the AGN in NGG5005. Despite this, the derived EW value places NGC5005 in 
the same regime as other Sbc morphological type galaxies in the Mb-EW and (Vmax)-EW planes 


®The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Insti¬ 
tute of Technology, under contract with the National Aeronautics and Space Administration. 
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as found in lLee et alJ (|2007l l. 


In addition to integrated properties from the aperture photometry, radial profiles were also 
created to examine the surface brightness values in each elliptical annulus. Figure 0] shows the 
radial changes in the surface colors B-3.6, R-3.6 ( 6 ), and B-R (c) (see also Section [2~ill . as well as 
the radial trend in EW (d). The radial trend in B-R is typical of large spiral galaxies, going from 
redder in the bulge-domina ted nuclear region to b luer in the disk where more of the current star 
formation is occurring fe.g.. iMarkarian et al.lll965l h After the initial gradient from red to blue, the 
B-R color remains roughly constant through the disk of NGC5005. There does not appear to be 
any change in B-R color corresponding to either the extent of the HI gas, R 25 , or the change in 
the outer disk exponential slope (Rbreak)- The radial trend in EW is much less regular than the 
B-R color trend. The higher values of EW correspond to spiral arms where there are more HII 
regions from recent star formation activity. Table [4] gives estimates of the EW gradient and B-R 
color gradient from a linear fit to these profiles. 

For visual reference, the characteristic radii described above are overplotted on both the Ha 
and R-band images in Figure [5j The images are displayed with the same scale to exemplify the 
concentration of the ionized gas emission (left) and extent of the stellar disk (right). As is expected, 
the morphology of the ionized gas emission in the central '^40" di rectly matches the morphology of 
molecular gas emission mapped in CO bv ISakamoto et al.l ( 1200 ( 1 ). As can be seen in the zoom-in 
image, this morphology features a central nuclear disk (first solid e llipse) with a ring of e mission at 
about 30^^-40^^ (second solid ellipse; referred to as the 3 kpc ring in ISakamoto et al.ll200nl ). The two 
central ellip ses have semi-major a xes of 10 ^^ and 40^^, respectively, assuming an inclination angle 
of 53° (from ISakamoto et al.ll2000l i and position angles of 75°. The Ha emission continues beyond 
this central region, tracing the spiral arms which appear to emanate from either end of the ring 
along the semi-major axis. The ionized gas emission is contained within the HI disk (first dotted 
ellipse). The R-band image on the right shows a high surface brightness disk component that 
extends beyond R 25 (middle dashed-dotted ellipse) followed by a low surface brightness component 
that begins slightly within the break radius (outer-most dashed ellipse). The HI disk (inner-most 
dotted ellipse) lies well inside the high surface brightness disk component. 


3. GAS KINEMATICS 

In this section, we return to the HI synthesis observations to explore the neutral gas kinematics. 
In addition, we use integral field unit (IFU) spectroscopy to extract ionized gas kinematics in the 
central region of NGC5005. We use both ionized and neutral gas kinematics to derive the circular 
rotational velocities of the gas as it follows the overall gravitational potential of the galaxy. Starting 
with the velocity field, we follow well established techniques to determine the bulk motions of the 
gas from which we extracted a rotation curve. This process is described in more detail below. 
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Fig. 5.— Left. Narrowband Ha image of NGC5005 from HDI. Right Masked and interpolated 
R-band image of NGC5005 from HDI. The inner two solid ellipses demonstrate the size of the 
nuclear disk and ring structure, respectively (see Section [2.5p . The outer three ellipses demarcate 
the three radii shown in Figure IH Rhi (dotted), R 25 (dashed-dotted), and Rbreak (dashed), listed in 
order of increasing size. 
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3.1. Neutral Gas Kinematics 


Beginning wi th the HI. we used stand ard analysis tools distributed as part of the GIPSY 
software package (jvan der Hulst et alJll992l l for analysis of HI synthesis data cubes to extract an 
intensity-weighted mean velocity field for NGC5005. To create this first moment map, individual 
channels from the low resolution datacube at 5 km s“^ channel”^ were smoothed by a factor of 3, 
and clipped at the 3 ct level, then interactively blotted to identify signal. Figure [T] (d) shows the 
hnal velocity held with isovelocity contours. The straight isovelocity contours seen in Figure [T](d) 
indicate that the rotational velocity of the gas remains fairly constant with radius with no obvious 
signs of a warp at large radii. Kinematic warps due to changes in position and/or inclination angle 
of the disk often result i n a v elocity held with a characteristic S-shaped distortion seen at larger 
radii (e.g., M83; lBosmalll98lh . The low signal-to-noise of the HI observations within the central 
~20" (1.6 kpc at 16.5 Mpc) makes it difficult to trace the HI gas motions. However, as introduced 
previously, I Sakamoto et al.l ((20001) found an S-shaped distortion in the GO velocity held. This 
central kinematic distortion is most likely due to the presence of a bar causing the gas to be in 
non-circular orbits. 

The hlled circles in the bottom panel of Figure[7]show the rotation curve for NGG5005 derived 
by htting a series of concentric tilted rings to the HI velocity held (GIPSY task rotcur). In this 
analysis the galaxy is described using a set of concentric rings, each with its own center coordinates 
(xq, Vo), systemic velocity (I4ys), expansion velocity, inclination angle (i), position angle (P.A.), and 
rotation velocity (Trot)- The radii and widths of the rings are dehned such that the rotation curve 
is sampled at a rate of two points per synthesized beam width, which in this case is every lO'h The 
expansion velocity is hxed to zero. Initial estimates for the center coordinates, inclination angle, 
and position angle were hrst obtained by htting ellipses to the HI intensity distribution out to 3a 
(GIPSY task ellfit). The systemic velocity was measured from the spectral prohle derived from 
the integrated intensity HI map. These initial values were fed into ROTCUR and iteratively ht to 
the receding and approaching sides together. Once the center position (xq, yo) was hxed, ROTCUR 
was run again to determine I^ys with i and P.A. as free parameters. Additional runs of ROTCUR 
were used with (xq, yo) and Fsys hxed to determine i and P.A. for a combination of both sides. 
The rings are allowed to vary in position and inclination angle as a function of radius to account 
for warps {top panels). However, point-to-point variations in i can be caused by effects such as 
st reaming motions alo ng spiral arms for which the htting routine tries to compensate by changing 
i (Ide Blok et al.l 120081 ). A smoothly varying or constant distribution in i and P.A. was adopted to 
derive the underlying bulk rotation and prevent spurious second-order effects in the rotation curve. 
Because variations in position and inclination angle were small, constant position and inclination 
angles were used for this analysis (dashed lines). Smoothly varying hts to both the inclination and 
position angle were also tried and determined to be unnecessary as this only changed the derived 
circular rotational velocities by a few km s“^. A hnal run of ROTCUR was done with all parameters 
hxed to derive the hnal azimuthally averaged rotation curve (bottom panel, filled circles). 


Uncertainties on the rotational velocities were conservatively estimated to be a combination 
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of non-circular thermal gas motions and inclination errors. The uncertainties on the rotational 
velocities for the first two radii are dominated by thermal gas motions and beam smearing. A 
velocity dispersion, or second moment map, was created from the high resolution HI data cube. 
Uncertainties due to thermal gas motions were estimated to be about 10% of the average velocity 
dispersion values, which is the accuracy to which we believe the line centroid can be estimated. The 
velocity dispersion for NGC5005 in the radial range 20'' to 40" is quite large; average values range 
between 40-50 km s“^. This is why the first two rotation curve points have noticeably larger error 
bars than the rotational velocities farther out where the dispersion drops to more typical values on 
the order of 15 km s“^. Outside of r > 40", the uncertainties are dominated by inclination angle 
errors. Due to the relatively high inclination angle (~68°) of NGC5005, an inclination angle error 
of ±10° in the inner parts and ±5° in the outer parts produces corresponding rotational velocity 
errors of <10% and <5%, respectively. Circular rotational velocities and their errors in the central 
20" were not derived from the neutral gas data due to the low signal-to-noise in the HI data. 

The rotation curve for NGC5005 is remarkably flat; it varies by less than 10 km s“^ over a 
range of radii from 20" to 140" (1.6 kpc to 11.2 kpc at 16.5 Mpc). The last measured rotational 
velocity is 265.2±22.0 km s“^ which yields a dynamical mass of (1.83±0.30) x 10^^ Mq at 140" (11.2 
kpc at 16.5 Mpc). If we assume that the rotation curve remains flat and extrapolate out to the 
break radius, this results in a dynamical mass of (2.82±0.47) x 10^^ Mq measured at 216.3" (17.3 
kpc at 16.5 Mpc). 


The HI rotation curve for NG C5005 was additionally deriv ed using a code for modeling asym¬ 
metries in disk galaxies (DiskFit: ISpekkens fc SellwoodI 120071 ). DiskFit uses a physically moti¬ 
vated model rather than a parameterization of concentric rings. It is capable of fitting models with 
both rotation and lopsided or bisymmetric no n-circular motions, as would be expected due to the 
presence of a bar (|Kuzio de Narav et al.ll2012l b Rotation curves for NGC5005 were derived using 
three different models: rotation only, lopsided flow, and bisymmetric flow. Due to the low signal- 
to-noise ratio of the HI at radii smaller than 20", only the rotation dominated parts of the velocity 
field were fit. Rotation curves derived from all three DiskFit models, therefore, were in good 
agreement with a rotation only model and the tilted ring fits from rotcur. Derived best fitting 
kinematic parameters from GIPSY’s ROTCUR and DiskFit’s rotation only model are presented in 
Table [5l 


3.2. Ionized Gas Kinematics 

Given the low signal-to-noise ratio of the HI observations in the inner 20" of NGG5005, we 
turn to ionized gas to provide kinematic constraints in this inner region. Figure [6] shows the ionized 
gas emission map as traced by the [NH]A6584 line flux (left) and luminosity weighted mean velocity 
field determined from cross-correlation of the entire spectrum (right). We used the [NHJA6584 
emission line to demonstrate the distribution of ionized gas in the central 70" of NGG5005, as the 
Ha flux is often affected by underlying stellar absorption in the diffuse gas regions. Regions where 
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Fig. 6.— Left: SparsePak [Nil] A6584 emission line flux map for the central 70" x 70" of NGC5005. 
The [Nil] A6584 flux map illustrates where the ionized gas emission is strongest and thus where 
the cross-correlation analysis has the highest signal-to-noise ratio. Right: Luminosity weighted 
mean velocity field for the central 70" x 70" from SparsePak observations of ionized gas. Contours 
represent lines of constant rotational velocity, the values of which are indicated above and on the 
left in units of km s“^. 


the [NII]A6584 line flux is strongest represent areas with higher signal-to-noise ratios with typical 
TDR values around 40 from the cross-correlation results. The lower signal-to-noise ratio areas still 
have reliable cross-correlation results with TDR values closer to 20. We used ROTCUR in GIPSY to 
derive a rotation curve from the ionized gas velocity field following the same steps for the HI data 
outlined above. The rotation curve derived from SparsePak ionized gas observations is shown in 
Figure [7] (open stars). The ionized gas velocity field was sampled every 5" out to a radius of 35". We 
allowed the position angle, inclination angle, and the systemic velocity to be fit by ROTCUR and not 
constrained to the values used for the HI derived rotation curve. Constant values of position and 
inclination angle were adopted (dash-dot lines). The ionized gas data fills in the central kinematics 
that are unconstrained by the HI observations and extends up towards the HI derived rotation 
velocity at 40". The inner two HI rotation curve data points display higher rotational velocities 
than the ionized gas, but this is most likely due to the effects of beam smearing. 


W hile the ionized gas morphology agrees well with the CO morphology found bv ISakamoto et al. 
(j200d L we do not observe as large of a velocity in the nuclear disk (r < 10"). This difference could 
arise because the ionized gas rotational velocities were derived from the centroid of the emission 
line detected in a 5" fiber and, therefore, do not include the high velocity dispersion of gas associ¬ 
ated with the AGN. We do, how ever, detect the same c entral S-shaped distortion in the ionized gas 
velocity field, as was detected bv ISakamoto et al.l (j2000l l in the CO. Given the similar morphologies, 
it is plausible that the ionized gas is subjected to the same bar-driven dynamics as the molecular 
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gas discussed extensively in ISakamoto et alJ (|200Cll i . 


4. DARK MATTER MASS DECOMPOSITION 


Rotation cu rve decomposition analysis is a powerful tool to probe the dark matter content 


in galaxies (e.g.. Ide Blok et al.l 120081 ). This well established analysis method combines the derived 


galaxy kinematics from the observed rotation curve with estimates of the baryonic mass distribu- 
tion in order to derive a n estimate of the dar k matter distribution as a function of radius (e.g., 
de Blok fc Bosmal 120021 ) . Ide Blok et al.l (120081 ) found that both widely used dark matter models 
(NEW and isothermal) fit equally well for galaxies with < -19. In this analysis, we fit only 
the spherical pseudo-isothermal dark matter halo model to the observations as the data are not 
sensitive to the central dark matter distribution. 

One of the strengths of the present analysis is the use of NIR data to derive stellar bulge 
and disk distributions. In particular, the Spitzer 3.6//m images are relatively insensitive to radial 
changes in mass-to-light (M/L) ratio that might be introduced by either internal dust obscuration or 
changes in the dominant stella r population. Indeed, the M/L r atio at 3.6/zm is nearly independent 


of star formation history fe.g.. lLi et al.ll2007l . iMeidt et al.ll2014l ). The combination of NIR, HI, and 


CO observations allow many of the uncertainties associated with the spatial distribution of both 
major baryonic components (gas and stars) to be minimized. 

The contributions from gas and stars to the overall rotation curve for NGC5005 were calculated 
from the gas mass and stellar light surface density profiles. The gas mass surface density was 
measured from the HI and CO integrated intensity maps. The gas distributions were derived 
separately and then added together to create a total gas mass surface density distribution. The 
contribution to the total circular rotational velocity was calculated from the to tal gas distrib ution 
using the GIPSY task rotmod assuming an infinitely thin gas disk. Following [Broeili (1 19921 1. the 
total gas mass was multiplied by a factor 1.3 to account for primordial helium. Figure [8] shows 
the gas mass surface density distributions for the HI and H 2 (as traced by CO) separately {filled 
circles and stars, respectively), as well as the resulting total gas distribution when added together 
{dashed line). The molecular gas is the dominant contribution to the total gas mass surface density 


Table 5. Kinematic Properties of NGC5005 


Parameter 

ROTCUR Value 

DiskFit Value 

Kinematic center (RA, Dec J2000) 

13:10:56.28, +37:03:31.7 

13:10:39.05, +37:05:35.4 

Systemic velocity 

948 km s“^ 

944 km s“^ 

Adopted position angle 

65.6° 

65.2° 

Adopted inclination angle 

68.1° 

67.8° 
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Fig. 7.— Rotation curves for VLA HI and SparsePak (SPak) ionized gas derived from tilted rings 
using the GIPSY task rotcur. The radial range ht for the HI was from 20"-l40" (1.6-11.2 kpc at 
16.5 Mpc) and 0 ^h 35 " (0-2.8 kpc at 16.5 Mpc) for the ionized gas. Top: Inclination angle fits to the 
approaching {filled circles) and receding {open squares) sides and the adopted constant inclination 
angle fit to both sides {dashed line) of the HI velocity field. Both sides of the SparsePak velocity 
field were fit at the same time {open stars), and a constant inclination angle was also adopted 
{dot-dash line). Middle: Position angle fits to the approaching {filled circles) and receding {open 
squares) sides and the adopted constant position angle fit to both sides {dashed line) of the HI 
velocity field. The open stars show the SparsePak position angles fit to both sides of the velocity 
field and dash-dot line shows the adopted constant position angle. Bottom: Circular rotational 
velocities fit to both sides for the HI {filled circles) and SparsePak ionized gas {open stars) derived 
using the adopted inclination and position angles at each radius. 
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for radii within 60'^ (4.8 kpc at 16.5 Mpc). The bump in the molecular gas distribution between 
30"-40" corresponds to the ring of emission observed in the narrowband Ha image (Section 12.51 
and Figure [5|). Model rotation curves for the stellar distributions were derived separately for the 
bulge and disk based on the decomposed light surface density profiles discussed in Section 12.41 The 
contribution to the total circular rotational velocity for these stellar components were calculated 
with ROTMOD. A vertical sech-squar ed distribution w i th a v ertical scale height of 0.93 kpc estimated 
using the fiducial relation derived in iBershadv et al.l (j2010|) was used for the stellar disk model. A 
spherical distribution was assumed for the bulge model rotational velocities. 

Figure [9] presents the rotation curve decomposition for NGC5005. Circular rotation velocities 
for the total gas, stellar bulge, and stellar disk were summed in quadrature and subtracted from 
the observed rotation curve. A spherical pseudo-isothermal dark matter halo model was fit to 
the residuals. The bulge and disk M/L ratios were left as free parameters to be fit during the 
decomposition. They were assumed to be constant with radius. For this galaxy, leaving the stellar 
disk and bulge M/L ratios in addition to the characteristic radius (Rc) and velocity (Vh) of the 
dark matter halo model as free parameters creates a poorly constrained problem without a unique 
solution. Therefore, the best fitting Rc, Vh, M/Lbuige, and M/Ldisk were achieved using interactive 
fitting. The best fit achieved in this ar ialysis resulted in a disk M/Ldisk = 0.5 and a bulge M/Lbuige = 
0.4. This disk M/L is consistent with Eskew et al. I (l2012l l who suggest that the M/L ratio at 3.6/rm 
should be about 0.5 based on stellar population models and measurements of the Large Magellanic 
Cloud. The bulge and disk M/L ratios result in a best-fitting dark matter halo model with a 
characteristic radius of 2.5±0.1 kpc and a velocity of 228±1.0 km s“^. The observed rotation curve 
could not be adequately described at large radii using the baryonic contributions alone (crosses in 
Figure [9|). 

The dip in velocity in the observed ionized gas rotation curve between 12"-24" (0.96- 1.92 
kpc at 16.5 Mpc) is an interesting feature which demonstrates the connection between the galaxy 
morphology and its kinematics. It appears to occur in the transition region between the nuclear 
disk and ring structu re, as seen in the Ha image (Figure [5] left) and detected in CO emission 
(jSakamoto et al.l 120001 ). This region corresponds to a decrease in the total gas contribution to the 
overall rotation curve, which is dominated by the molecular gas at these inner radii (Figured]), as 
well as to where the stellar contribution to the overall rotation curve transitions from being bulge- 
to disk-dominated. It is also visible as a prominent dip in the equivalent width profile seen in 
Figured] (d). 

The mass decomposition reveals that NCC5005 is baryon dominated at the inner-most radii. 
NCC5005 follows the same trend as most other massive, luminous galaxies where its mass is baryon 
domi nated in the core an d gradually becomes more dark matter dominated with inc reasing radius 


(e.g.. Ide Blok et al.M2n08l ). Furthermore, galaxies of similar dynamical mass from Ide Blok et al 


(|2n08l L such as NCC6946 and NCC7331, display rotation curve decompositions similar to that of 
NCC5005 where the dark matter halo contribution to the total rotation curve does not overtake the 
stellar disk contribution to the total rotation curve until close to or beyond the last measured point. 
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Fig. 8.— Atomic and molecular gas mass surface density distributions for NGC5005. The atomic 
(filled circles) and molecular (stars) gas distributions are derived from the HI and CO integrated 
intensity maps, respectively (see Figure [U (a)). The spatial sampling for each distribution is roughly 
half the beam size from their respective maps (20" for HI and 6" for CO). The dashed line represents 
the total gas mass surface density when the atomic and molecular distributions are added together. 
Mass surface densities are plotted as the log of the value. 
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The fractional contribution of the maximum stellar disk to the rotation curve at 2.2 times the disk 
scale length is only 0.71 in NGC5005. This implies that the stellar disk can only explain about 
70% of the observed rotation curve at the point where the stellar disk’s contribut ion should be 


great est. This is slightly less than average values found for other samples of galaxies (|Swaters et al, 


201 ll and references therein), which were between 0.8 and 0.9. 


5. CONCLUSION 

We have acquired multiwavelength observations to examine the gas and stellar content in 
NGC5005. We have utilized new VLA G configuration observations to explore the neutral hydrogen 
gas morphology and kinematics, as well as SparsePak IFU observations to derive kinematics of the 
ionized gas in the central region of NGG5005. Deep 3.6/xm images from Spitzer allow us to trace 
the extended stellar population. Finally, we have used optical broadband B and R and narrowband 
Ha observations to detect properties of the dominant stellar population and star formation activity. 
The main results from this multiwavelength study are summarized below. 


1. The most striking morphological feature of NGC5005 is its extended stellar population. The 
neutral and molecular gas components are both conhned well within the high surface bright¬ 
ness stellar disk. The low surface brightness stellar component appears to extend nearly 3 
times beyond the gas disk detected in HL The central 40^^ of NGC5005 has a lack of detectable 
HI, but displays both molecular and ionized gas emission. 


2. Spectra of the ionized gas in the central region of NGC5005 reveal signatures of both thermal 
and non-thermal emission processes. The nuclear AGN region displays broad emission lines 
with non-thermal line intensity ratios. Double-peaked emissio n lines are seen to the northwest 


of the nucleus, approximately where iSakamoto et al.l (120001 ) found dynamical evidence for 


infalling gas potentially due to an inconspicuous stellar bar. The ionized gas emission outside 
of the nuclear disk appears to be dominated by thermal processes from recent star formation 
activity. 


3. Inspection of the optical broadband images reveals numerous known and suspected com¬ 
panions of NGC5005. These include SDSS J131051.05-1-365623.4 which was detected in the 
VLA observations and SDSS J131105.57-1-371036.1 which appears to be possibly connected 
to NGC5005 by a faint stellar stream. 


4. Surface photometry of NGC5005 shows a NIR surface brightness distribution with a central 
bulge component and an inner and outer exponential disk. The radial prohle of optical color 
displays a distribution that is typical for large spirals like NGC5005, and is consistent with 
NGC5005 having little current star formation activity. The unusually low star formation rate 
and equivalent width measurements from the Ha imaging are expected given the gas poor 
nature of NGC5005. 
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Fig. 9.— Decomposition of the rotation curve from the HI, optical, and NIR data. The total gas 
{filled squares), stellar bulge {open triangles), stellar disk {open stars), and dark matter halo model 
{dashed line) components are added in quadrature to achieve the best overall fit {solid line) to the 
observed rotation curve. Open circles represent circular rotation velocities derived from the ionized 
gas (SparsePak) observations and the filled circles are from the HI (VLA) observations. Adding 
the baryonic components in quadrature without the dark matter results in the crosses. This total 
baryonic contribution primarily follows the stellar distribution in this gas-deficient galaxy. 
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5. A rotation curve for NGC5005 was derived from both the neutral and ionized gas kinematics. 
The ionized gas rotation curve in the central region of NGC5005 exhibits structure which 
can be connected directly to the ionized gas morphology (which follows the molecular gas), 
and also appears to correspond to the stellar distribution as it transitions from bulge to disk 
dominated. The HI rotation curve at large radii is remarkably flat and exhibits no signs of 
warps. 

6. The gas kinematics and surface densities were used in conjunction with the stellar bulge and 
disk light profiles to decompose the observed rotation curve of NGC5005 and derive a model 
dark matter distribution. The results of the mass decomposition demonstrate that NGG5005 
is baryon dominated in its central region and that the prominent stellar disk only maximally 
contributes about 70% to the total rotation curve. 
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